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^ ' Abstract. We present a Green's function based treatment of the effects of electron-phonon 

coupling on transport through a molecular quantum dot in the quantum limit. Thereby we 

combine an incomplete variational Lang-Firsov approach with a perturbative calculation of the 

2 ' electron-phonon self energy in the framework of generalised Matsubara Green functions and 

a Landauer-type transport description. Calculating the ground-state energy, the dot single- 
^- particle spectral function and the linear conductance at finite carrier density, we study the low- 

'-rJ I temperature transport properties of the vibrating quantum dot sandwiched between metallic 

S ■ leads in the whole electron-phonon coupling strength regime. We discuss corrections to the 

f^ ' concept of an anti-adiabatic dot polaron and show how a deformable quantum dot can act as a 

molecular switch. 



ly-^ , 1. Introduction 

lO ' Recent progress in nanotechnology allows for the fabrication of electronic devices with 

l/^ . organic molecules as the active elements, which may constitute an alternative to conventional 

^P I semiconductor technology in the search for further miniaturisation. The basic example for 

such a device is a single organic molecule contacted with metallic leads. The molecule can 
be described as a quantum dot, i.e. as a system of finite size that is coupled to macroscopic 
charge reservoirs. The electronic transport properties of such molecular electronic components 
^ ■ depend on the geometry of the molecule and the properties of the molecule-lead contact, as well 

as on the average charge of the molecule [1]. Since the molecule is so small, quantisation of 
energy levels becomes important. Furthermore, the molecule is susceptible to specific structural 
changes in the presence of charge carriers. The energy of these deformations or vibrations 
can be comparable to the kinetic energy of passing particles, whose mobility may therefore be 
substantially affected. This process seems to play a fundamental role in the observed non-linear 
behaviour of basic molecular devices, including negative differential resistance and hysteresis [2]. 
With strong coupling to vibrational degrees of freedom, current rectification and fast current 
switching may be realised [3]. 



2. Model 

To investigate transport in a molecular quantum dot, we consider the model Hamiltonian 

^ = E ^fc4aCfca + ^dU - guoib^ + b)dU + u:ob% - 4= E {d^^ka + cLd) . (1) 

k,a k,a 

Here c^^ (c^a) are creation (destruction) operators of non-interacting electrons with energy 
Ek [k = 1, . . . ,N) in the left and right lead (a = /,r). We assume that the leads are semi- 
infinite one-dimensional chains with a semi-elliptical density of states q{^) = -^ ^^ 6{^ — Ek) = 
{2/TTW^)^yW'^ - ^2 Q{W'^ - ^2), where W is the half bandwidth. The quantum dot is modelled 
by a single level A with fermionic operators S^'. Working with spinless fermions, we take a 
large local Coulomb repulsion at the quantum dot for granted. To describe the deformation of 
the molecule, an electron at the dot couples via a Holstein-like term (ex g) to a local phonon 
mode b^^' of energy loq. The last term in ([T]) allows for the dot-lead particle transfer (ex 1^). 

In analogy to Holstein's small polaron theory [3], we expect - for sufficiently large electron- 
phonon (EP) coupling g and frequency uq ~ the formation of a polaron-like state at the dot, 
corresponding to an electron with an accompanying phonon cloud. To account for this effect, a 
generalised Lang-Firsov-transformation H = UHU^ with U = ex.p{g{b^ — b)d^d} and g = ^g is 
applied, where 7 G [0, 1] is a variational parameter. The transformed Hamiltonian reads 

H = ^E,clc,^ + AdU -C,dU + uobh-Y,{Ctd^c,^ + C}cld) , (2) 

k,a k,a 

with the renormalised dot-level A = A — ep7(2 — 7) and the polaron binding energy Sp = g'^ojQ. 
As 7 grows from zero to one, our approach interpolates between a weak coupling ansatz (7 = 0, 
H = H), and the complete Lang-Firsov-transformation (7 = 1) restricted to large phonon 
frequencies and strong EP coupling, when the direct EP coupling via C^ = 5Wo(l — 7)(^^ + b) is 
replaced by a phonon-affected dot-lead transfer term Ct = {td/y/N) exp{— ^(6^ — 6)}. In this way, 
we are able to describe the system for a large range of parameter values. Note that the effects 
of a finite Coulomb interaction might be included by considering instead of the Hamiltonian p!|) 
Hubbard/Anderson-Holstein-type models O [6] . Here the influence of the phonons is mainly to 
suppress the repulsion between the electrons at the molecular orbitals [7]. When the energy 
scales set by the Coulomb and EP interactions become comparable bipolaronic states may form 
at the quantum dot [5] . 

3. Theoretical approach 

Our main interest lies in the single particle spectrum at the dot and the linear conductance in the 
case of vanishing voltage bias between the leads. Both of these quantities can be obtained from 
the polaronic spectral function ^^^(a;) = lim5^o+[ImG^^(ti; + i(^)], where the retarded Green 
function G^ corresponds to the (polaronic) operators d^"^' in the transformed Hamiltonian in 
equation ([5]). We base our calculations on the equations of motion of generalised Matsubara 
Green functions in equilibrium [51 E] 

Grfd(ri,r2;{y}) = -^{Trd{n)d\r2)S) . (3) 

The mean value and the time dependences in Q are determined hy H — fiN, where fi is the 
equilibrium chemical potential of the system and N denotes the particle number operator. The 
S-matrix 

5 = r.exp|-j^ dry, (T)Q(r) + V;(T)C/(r) + y,(r)C,(T)| (4) 

describes the coupling to the components of a fictitious external potential {V} . 



3.1. Dot spectral function 

Starting from the equations of motion for the dot Green function Gdd and the lead-dot transfer 
Green function Gcd\ka, we express the occuring interaction terms by functional derivatives with 
respect to {V} [8], e.g., 

^^^^ -Ctin,{V})G,d-Min,r2,{V}) + ^^^^ , (5) 

with the interaction coefficients C^ {t;{V}) = (T^C^ (r)5)/(5). The resulting coupled 
equations lead to afunctional differential equation for the polaronic self energy S^^(ri,r2; {V}) = 
G^J {ti,T2) — G~^^ ("T"!, 7"2; {y}) (see equation (17) in [9]), which can be evaluated in a two-step 
process. First, we neglect terms with functional derivatives of S^^. We then insert the result, 
EJ^^ , into these derivatives and only keep terms up to second order in C^ . In this way, we find 

42(n,r2;{F}) = -C,(ri;{FMri-r2] + VA(n;{n)GSfc,(ri,r2)Cl(r2;{F}) (6) 



EGS.a(n,r2) 
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_{S) 
+ G)IJ{n,T2) [{TrCd{n)Gd{T2)) - Cain; {V})Cd{T2; {V})] . 



Here G^J denotes the Green function determined by the ffi'st order self energy and G^Jf^^ is the 
free Green function of the lead states. We then let {V} — )■ and calculate the correlation 
functions of the interaction coefficients assuming an independent Einstein oscillator. After 
Fourier transformation and summation over bosonic Matsubara frequencies, the retarded self 
energy follows in the low-temperature approximation f3ujQ 3> 1 as [9] 

/w -I 
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Within our iterative scheme, the corresponding spectral function A^J = Add is evaluated in a 
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where the first order self energy T,^J is given by the first two terms in ([7]). The self energy 

(2) 

S^^ = S^^ accounts for multi-phonon processes as well as finite particle densities. As 5 — )• 0, 
our model ([1]) reduces to the Fano-Anderson Hamiltonian for a rigid impurity in a ID lattice. 
Then only the first term in ([7]) remains and gives the exact self energy. For finite EP coupling, 
the spectrum contains multiple phononic side bands. If the condition ujq < W — \n\ is fulfilled, 
these bands overlap and lmT,dd 7^ along the whole cj-axis. Otherwise, a localised polaron-like 
state may exist in the intervals where ImS,^;^ = 0. 



For the numerical evaluation of the spectral function we keep (^ in (l8|) as a small positive 
parameter (5 < 5 • 10~^). This avoids the problematic evaluation of principal value integrals in 
HeTidd for (5 — )• 0. The resulting spectral function fulfils the sum rule J Add{'^)d^ = 1, which is 
preserved in our approximations. 

3.2. Ground state energy 

The variational parameter 7 is determined by minimisation of the ground state energy E = 
[H — fiN) with respect to 7. We identify the statistical averages in E with expressions like 
equation ([5|), letting t^ — )• Tg" and {V} — )■ 0. As in the evaluation of S)^^ , we neglect the 
functional derivatives in equation ([5]), so that {C^Scf,^ k, —{C^)Gcd;ka{'Ti,T2)\^ -i-r^' ^PP^yi'^S 
the same approximation to the equations of motion of the Green functions, we determine their 
Fourier transforms in the complex plane to first order as 

G..,Uz) - Gf}^,Sz) + {Cl){C,)[G^}^,^{z)]'Gdd{z), (9) 

GcdMiz) « -(C/)G(2,,(z)Grfrf(z) . (10) 

For the variation of 7 we can omit the constant lead energy related to G^^ fca- -'■^ equilibrium, 
(6'6) and {Gd) vanish as T — t- 0, and we obtain 



E = -2tle-3 / diQii) / doo'Add{oj') 
J-w Jo 

rW rO 

+ 2tle-3 / d^^(e) / div' Addioj') 



+ 



-(e-/^) (^'-(e-/^))2 
2 ^ i-ii 



(e-/x) (^'-(e-M))^ 



^ J —00 

+ (K-p)[ dJAddioj'). (11) 

J —00 

Via Add, ^ and 5, the energy E, given by equation dH}, depends on 7 and can be used to 
determine the extremal variational parameter 7mm- 

3.3. Conductance 

We finally calculate the linear conductance L starting from the Meir-Wingreen formula |10j . 
which expresses the current J through the dot in terms of the retarded non-equilibrium Green 
function G^ of the electronic operators d in ([T|). For equal coupling to the left and right lead, 
and finite voltage bias $ = — (/x; — /ir-)/e, the current takes the simple form 

J=-etl / diQ{i)[nF{i-^Xl)-nF{i-^ir)]l^GUi)■ (12) 

J-w 

Here np{^) = (e^^ + 1)~^ is the Fermi distribution function of the isolated leads in thermal 
equilibrium, at respective chemical potential //;/,,. For small voltage bias, i.e. ^i^r = /^ i (^/i/2, 
we can express the current as J = —L6^i/e. Then, the linear conductance L = lim5^_!,o{— eJ/5/i} 
results from (fT2]) as 

j-w 

L = e\tl / de£'(e)[-n'^(e-A^)]^^,(e-/^), (13) 

J-w 

where the electronic spectral function A'^^{uj) = lim5_j.o+[ImG^^(a; + 16)] is now calculated 
in equilibrium. Based on the factorisation of the statistical averages of phonon and polaron 
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Figure 1. For model parameters ujq = 0.05, A = 1, t^^ = 0.25, // = 0. Panel (a): Linear 
conductance L and occupation of the dot n (inset) as functions of £p. At £p = 1.03 an abrupt 
transition takes place. Panel (b): Dot energy £■ as a function of 7 for several EP coupling 
strengths. At eL a second minimum at 7 = 1 becomes the global minimum, causing a jump in 
the extremal parameter 7min. (inset). Panels (c) and (d): Spectral function Addioj), integrated 
spectral weight <S'(i^) = J_^ du' Add{i^') and imaginary part of the self energy T,dd{^) for Sp in 
the vicinity of the transition, showing the sudden formation of a long-living polaron-like state. 

variables, a relation between A^^ and the polaronic Add has been derived in equation (40) of 
^ddi^) = e"^ Yl i~ [^'^'^('^ ~ ^^o)Q{^ - sujo) + AddiuJ + sujo)Q{-uj - sujq)] . (14) 

s>0 

Then, from (J13p and ()14p . the linear conductance in the low-temperature approximation is 

L = e\tle-^'gifi)AddiO). (15) 

As usual, it depends on the accessibility of dot states at the Fermi level. In addition the 
conductance exhibits a Lang-Firsov renormalisation td = tde~^ '"^ of the dot-lead transfer 
integral, in accordance with the results for the Holstein model. 

4. Numerical results and discussion 

With our choice of a semi-elliptical lead density of states, the non-interacting system (cp = 0) 
is translational invariant if the dot-lead transfer td is equal to the value t = W/2 for electron 
hopping along the leads. In the following, we fix t = 1 and set td = 0.25, representing the weak 
dot-lead coupling case. 



4-1. Adiabatic regime 

Let us first consider the adiabatic regime (loq = 0.05 < td), where the phononic time-scale is 

much slower than the electronic time-scale and the deformation of the dot adjusts quasi-statically 
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Figure 2. For model parameters ujq = 3, A = 1, td = 0.25, /x = 0. Panel (a): Because E has a 
single minimum at 7mm for all £p, a smooth crossover with 7mm > 0.7 occurs (inset). The linear 
conductance reaches a maximum at Ep ~ A, when the peak in Add is shifted to the Fermi level 
(c.f. panel (b)). Panel (c): For large £p, a bound polaron-like state may exist in the intervals 
where ImS^d = 0. Then the electronic spectrum exhibits a multi-peak structure. 



to the average electronic occupation. We set /i = and A = 1, so that without EP coupling 
(ep = 0) the dot simply acts as an impurity leading to scattering. For these parameters the 
conductance L = 3 • 10~^ is small, i.e. far below its value L ~ 0.16 for A = (see figured^). As 
Ep increases, L and the dot occupation n = (Sd) grow only marginally. At a critical EP coupling 
Ep = 1.03 a sudden transition to n « 1 takes place, resulting in a drop of L by eight orders of 
magnitude. This behaviour can be traced back to a jump in 7mmi as figure [Tja shows. For small 



Ep, E{^) has a single minimum at 7,1 



0.02 ^ 1. Hence there is almost no renormalisation of 



A or td and L is nearly independent of Ep. However, the third term in ([7]) contributes and the 

single peak in A^^J is added to ImS^^, resulting in two maxima of equal spectral weight around 
w = A in Add (see figure [It). At Ep = Ep, when the gain in potential energy overcompensates 
the loss in kinetic energy, a second minimum of £^(7) appears and becomes the global minimum. 
As a result 7mi?i, jumps to unity, which corresponds to a complete Lang-Firsov transformation. 
The dot level (A = A — e^ = —0.03) is shifted below the Fermi level and the dot-lead transfer 
is reduced by a factor of exp{— ^} with Ip > 20, so that L vanishes. Because ujq <^ W, 
the spectrum consists of overlapping phononic bands around a pronounced central peak, which 
signals the formation of a long-living polaron-like state at the quantum dot. We know from 
variational approaches to the polaron problem, that the jump in 7mm niay be an artefact of 
our variational ansatz. In the Holstein model with EP interaction at every lattice site, the 
formation of a heavy polaron can occur as a sharp but always continuous crossover for small 
phonon frequencies. In contrast, a true phase transition from zero to finite dot occupation n is 
found for a single electron at the quantum dot with A > |12j . This phase transition becomes 
more pronounced as ojq gets smaller. In [5], similar behaviour is found in a generalised Holstein- 
Hubbard model with site-dependent potentials and EP couplings. Therefore, our results do not 
contradict the overall physical picture, and the variational Lang-Firsov transformation simulates 
the rapid adiabatic transition by a discontinuous change in 7mm- 



4.2. Anti- adiabatic regime 

In the anti-adiabatic regime with cjq = 3 ^ td, the dot deformation adjusts instantaneously to 
the presence of an electron. As figure [2^ shows, our approach is better suited to this regime: 
For /.i = and A = 1, E('y) has a single minimum 7mm > 0.7 for all EP coupling strengths and 
we find a smooth crossover to n ~ 1. At the quantum dot, both the renormalisation of td and 
A affect the transport. Since the phonon number g^ = Ep/ujQ is small in the anti-adiabatic case, 
the shift of the dot level is the predominant effect. With increasing Ep, the repulsive potential 
A is continuously lowered and the conductance L grows, reaching almost the maximum of the 
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Figure 3. For model parameters cjq = 0.2, td = 0.25, /.i = 0. Panels (a) and (b): L and 'jmin 
as a function of the quantum dot level. For small Ep = 0.2, we find a smooth crossover. As Ep 
grows, a jump-like transition occurs and the dot acts as a molecular switch. Panel (c): Spectral 
function at the point of maximum conductance in the smooth crossover. Panel (d): Due to 
overlapping phononic bands, no bound state exists even for large Ep. 

A = case when Ep w A. Here the peak in ^4^^ is shifted to the Fermi level (see figure [2)3) and 
we can speak of phonon-assisted transport. Because ujq > W, the self energy features few non- 
overlapping phonon bands (see inset of figure[2]3). For large Ep, a bound polaron-like state forms 
when A is located in the intervals with Im Ti^d = 0, as figure[2}: shows for Ep = 3.5 and A ~ —2.5. 
The electronic spectrum then exhibits the typical structure of multiple Poisson-weighted peaks 
(inset). Note that the delta peaks have finite width due to our numerical parameter 6 > 0. 

The results show that, as for the Holstein model, the complete Lang-Firsov transformation 
is restricted to large phonon frequency and strong electron-phonon coupling. Away from this 
limit, 7mm significantly differs from one, and our approach allows for important corrections in 
the regimes of weak coupling and moderate to small phonon frequencies. 



4-3. Intermediate phonon frequency 

Next we investigate the intermediate regime, where phononic and electronic energies become 
comparable (cjq = 0.2). We consider an experimentally relevant situation by keeping the EP 
coupling fixed and varying the dot level (gate voltage). In the case of a rigid quantum dot 
{Ep = 0), scattering off the dot potential reduces the conductance. Therefore, in figure [3^, L 
exhibits a maximum at A = and decreases symmetrically as |A| grows. For finite EP coupling 
Ep = 0.2, this maximum is shifted to A > 0, where the repulsive dot potential is compensated 
by the EP interaction. The spectral function then shows a peak at the Fermi level and two 
phonon satellites (see figure [3J3). As we see from figure [3}3, the optimal variational parameter is 
a continuous function of A taking values from 0.1 to 0.9. Therefore, the effective renormalisation 
of td and A depends, via 'jmin, on the dot level A itself. In contrast to the result for a complete 
Lang-Firsov transformation with fixed 7=1 (cf. figure 5 of [13]), the conductance maximum 



is shifted by less than £p, and L decreases asymmetrically away from this point. In accordance 
with Galperin et al [13] and Mitra et al jl4) . we find no phonon side peak in L{/\). If £p is set 
to larger values {sp = 0.4 and 0.5), a transition in the dot occupation appears which is again 
related to a jump in 7mm- La Magna and Deretzis |2] found a similar behaviour: for intermediate 
phonon frequencies and strong EP coupling their variational ansatz showed bistabilities causing 
a sudden occupation of the quantum dot (cf. figs. 2c and 2d in [2]). As the dot level falls below 
a critical value, the charging of the quantum dot is accompanied by a drop in the conductance. 
In this way the quantum dot acts as a simple molecular switch. Now the spectrum shows a 
pronounced maximum at A, which is below the Fermi level. Note that, due to overlapping 
phonon bands, no bound state exists even for large £p (figure [H]l). 



5. Conclusion 

In the present contribution we studied transport through a deformable quantum dot, realised 
e.g. by an organic molecule in contact to leads. Our treatment is based on Green's functions 
equations, with the dot Green's function as the central quantity of our calculation. EP 
interaction significantly affects the shape of the associated dot spectral function. Since at low 
temperature the conductance of the quantum dot is determined by the spectral weight close to 
the Fermi energy, it changes accordingly. The basic effect, with dramatic consequences for the 
conductance, is the formation of a "localised" polaron-like dot state. Similar to the Holstein 
polaron, this effect can be captured in a Lang-Firsov approach. The virtue of our analytical 
approach lies in the variational determination of the Lang-Firsov parameter. This allows us to 
account for the basic polaronic effects both in the adiabatic and anti-adiabatic regimes, which 
are distinguished by their different influence on the spectral function and, accordingly, on the 
conducteance (briefly summarised in figure H]). Particularly for comparable phonon and electron 
time-scales we find interesting physical behaviour, which motivates further investigation of the 
quantum dot system, e.g. with respect to applications as a current switch. 

The limitations of our study suggest two directions for improvements. First, the inclusion 
of Coulomb interaction at the dot. Many interesting effects we present here arise from the 
competition between electron-phonon interaction and a repulsive dot potential. Coulomb 
interaction gives rise to an effective repulsive dot potential, whose strength depends on the 
dot population. This may allow for even stronger non-linear behaviour, as then the transition 




Figure 4. Sketch of the basic polaronic effects on transport through a vibrating quantum dot. 
In addition to a renormalisation of the dot-lead transfer integral, in the anti-adiabatic regime, 
a continuous lowering of the dot potential towards the Fermi level enhances the transport. In 
the adiabatic regime, the sudden formation of a quasi-localised polaron state causes a drop of 
the conductance. 



between an effectively attractive or repulsive dot depends on the dot population through both 
interaction mechanisms. Second, all our approximations neglect correlations between phononic 
and polaronic degrees of freedom (note however that the Lang-Firsov transformation introduces 
strong correlations between the phononic and electronic degrees of freedom). In this way, the 
variational Lang-Firsov transformation can be used to map the entire problem approximately 
onto a purely electronic problem, as in reference |2]. This reference reports results similar to ours, 
e.g. for the dot population, but also for articifial bistabilities corresponding to different local 
minima of the groundstate energy, which are a common feature of these variational approaches. 
As a first step towards the reintroduction of electron-phonon correlations our treatment includes 
the weak-coupling contribution to the self-energy, which partly accounts for retardation of 
electron-phonon interaction at finite phonon frequencies. A detailed comparison of our results to 
those presented in reference [2] has to be given elsewhere, and the improvement of our treatment 
along the lines indicated is the subject of future work. 
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